Ultracold atoms in one-dimensional optical lattices 
approaching the Tonks- Girardeau regime 
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Recent experiments on ultracold atomic alkali gases in a one-dimensional optical lattice have demonstrated 
the transition from a gas of soft-core bosons to a Tonks-Girardeau gas in the hard-core limit, where one- 
dimensional bosons behave like fermions in many respects. We have studied the underlying many-body physics 
through numerical simulations which accommodate both the soft-core and hard-core limits in one single frame- 
work. We find that the Tonks-Girardeau gas is reached only at the strongest optical lattice potentials. Results for 
slightly higher densities, where the gas develops a Mott-like phase already at weaker optical lattice potentials, 
show that these Mott-like short range correlations do not enhance the convergence to the hard-core limit. 

PACS numbers: 05.30.Jp, 03.75.Hh, 03.75.Lm 



Introduction. — Since the prediction by Jaksch et al. |1] 
on the experiment by Greiner et al. ultracold atoms in 
optical lattices have been the focus of much activity. By 
tightly confining the motion in the transverse direction, an ar- 
ray of quasi one-dimensional optical lattices results Q], where 
particle exchange between the one-dimensional tubes is sup- 
pressed. The role of quantum fluctuations is enhanced in one 
dimension compared to the three-dimensional case, such that 
traditional mean-field theories fail. Instead, the long-range 
low-energy physics is described by the Luttinger liquid model. 
In the limit of infinite repulsion between the atoms the atomic 
gas is called a Tonks-Girardeau gas f^, Isl] (TG). Because of 
the blocking of double occupancies, the resulting hard-core 
bosons have some properties very similar to non-interacting 
fermions: e.g. the density profiles become indistinguishable. 
However, more complicated properties such as the momentum 
distribution remain discriminating characteristics |6]. The 
regime of strong repulsions between bosons has been studied 
experimentally |7] and theoretically |8] for an atomic gas not 
subject to an optical potential, but the acquired values for the 
ratio of the repulsive interaction strength to the kinetic energy 
were rather low and the TG regime was not seen. By using 
an optical lattice, much higher values for this ratio could be 
reached |6]. The interpretation of these experiments is com- 
plicated by the finite-size effects due to the harmonic trap. 
But even in the homogeneous case, an accurate theoretical de- 
scription of the transition from a weakly interacting Bose gas 
to a strongly interacting Tonks gas has to rely on numerical 
simulations. Our aim is to model the experimental results of 
Ref. 1 6] using one single numerical framework which accom- 
modates both the weakly and the strongly interacting regime. 

The physics of ultracold atoms in optical lattices can be 
described by the Bose-Hubbard model |1], which considers 
bosons occupying Wannier orbitals. The validity of this model 
is confirmed by the ratio of the central to first Bragg peak in 
the experimentally observed momentum distributions, which 
depends only on the shape of the Wannier orbitals (see below). 
The TG regime is characterized by the absence of double oc- 
cupancies in the many-boson wave function. To identify the 



TG regime unambiguously, one has to evaluate whether the 
experimental results are better described by soft-core bosons 
with a considerable overlap or by hard-core bosons for which 
double occupations are explicitly suppressed. Exact results 
for realistic parameters over the entire range of the axial op- 
tical lattice depths used in the experiment are obtained us- 
ing quantum Monte Carlo methods. We find that the results 
for soft- and hard-core bosons do not coincide except for the 
strongest optical potentials used in the experiments, in con- 
trast with the fermionization approach of Ref. |6] which as- 
sumes hard-core bosons at all optical-potential strengths. 

One-dimensional optical lattice. — When an ultracold Rb 
gas of atoms is cooled and loaded into an optical lattice 
with very tight transverse confinement, its dynamics is gov- 
erned by the one-dimensional Hamiltonian, 
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with m the atomic mass, x,- the position of atom /, Vo(x) = 
Vo sin^ the optical potential (Vo takes the laser intensity 
and the dynamic polarizability of the atoms into account) and 
Vt{x) the harmonic trapping potential, which varies slowly 
compared to the optical potential. The wave vector k of the 
laser along the axial direction defines the length scale A,/2 
through k — iTi/'k and the recoil energy m) 
which we will use as an energy scale. The Hamiltonian Eq.([0 
reduces to the exactly solvable Lieb-Liniger |5] model for 
Vt{x) = and Vq{x) = while it reduces to a Mathieu equa- 
tion for Vt{x) = and = 0. The interaction g„„ between 
the atoms is determined by the three dimensional scattering 
length as of the atoms. Olshanii |9] studied the scattering 
problem of two particles in tight waveguides and found for 
the effective one-dimensional coupling constant 
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where = y^h/m(0± is the characteristic length of the trans- 
verse harmonic confinement. For very tight radial confine- 
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ment it suffices to integrate over the y and z directions assum- 
ing harmonic confinement, yielding gmt = ' . The Wan- 

nier orbitals are calculated for the periodic potential given 
by the kinetic and the optical terms in Eq.([0, restricted to 
the lowest band 1 10]. For low-density gases we can express 
the Hamiltonian of Eq.lITl in the Wannier basis, resulting in a 
Bose-Hubbard model |1, 11], 

h^-jZ + y L «'(«<■ - 1), (3) 

where the first summation runs over nearest neighbors only, 
the operator n, counts the number of bosons at site / and 
the effective parameters f/,e, and J represent the strengths of 
the on-site repulsion, the harmonic trapping and the kinetic 
hopping, respectively. Recent studies of the one-dimensional 
Bose-Hubbard model mainly focused on the Mott-superfluid 
transition, using a wide range of methods: a slave-boson ap- 
proach I12II . the numerical renormalization group fisl . the 
density matrix renormalization group ] 14], the time-evolving 
block decimation method ] 15] and Monte Carlo methods ] 16]. 
The main uncertainties in the model relate to the accuracy of 
the scattering length and the renormalization of the effec- 
tive parameters of the Bose-Hubbard model. As we work in 
the grand-canonical ensemble, the chemical potential /j must 
be fine-tuned such that the expected number of particles cor- 
responds to the experimental number of particles. The Bose- 
Hubbard model is simulated using the stochastic series expan- 
sion method ] 17] (SSE) with locally optimized directed loop 
updates ] 18, 19]. From this Monte Carlo simulation thermo- 
dynamic observables such as the energy, the (local) density, 
the (local) compressibility and the one-body correlation func- 
tion can be computed exactly in a statistical sense ]20] . 

Homogeneous system. — First, we consider a homoge- 
neous (£, = 0) atomic gas in a lattice of L = 128 sites with 
periodic boundary conditions at a very low but finite tempera- 
ture, T /J ~ 0.2. In Fig.^we show the internal energy per site 
of this system for increasing values ofU/J, keeping the aver- 
age density fixed at («} w 0.5. An ideal Bose gas occurs in the 
limit of vanishing U, which is indicated by the lower horizon- 
tal line in Fig. [2 while the ideal Fermi gas is found for U 
and indicated by the upper horizontal line. For very large val- 
ues of U /J, no site of the lattice will be doubly occupied and 
one can apply the Jordan-Wigner transformation to map the 
bosons onto fermions 1I21I1 . For small but finite U {U j J 
in Fig.Q , the system is adequately described by the standard 
Bogoliubov approximation ]10]. For large U a perturbation 
of interacting fermions was derived in Ref. ]22] to order \/U . 
From the log scale in Fig.[2it appears that the limit of non- 
interacting fermions is reached slowly for values of t//7 > 10. 
For U /J ^ I one has to resort to numerical methods, and we 
see that the SSE method remains efficient over the entire U/J 
range. Higher temperatures lead qualitatively to the same re- 
sults, but the description in terms of fermions is only valid for 
higher values ofU/J. Temperature can be seen as a source for 
exciting double occupancy on a particular site, whose likeli- 
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FIG. 1 : Internal energy per site as a function of U /J for a homo- 
geneous model with 128 sites at a temperature T/J = 0.2. The data 
points with error bars connected by the full line are the energies ob- 
tained by the SSE method. The energies for non-interacting fermions 
(upper dotted horizontal line) and non-interacting bosons (lower dot- 
ted horizonal line) are shown, together with the Bogoliubov approx- 
imation for bosons (dashed line on the left) and a first order pertur- 
bation theory for fermions (dashed line on the right). 

ness must be suppressed by a stronger on-site repulsion term. 

Inhomogeneous system. — The harmonic trapping po- 
tential breaks the homogeneity of the system. For the pa- 
rameters we follow Ref. f^: the scattering length of Rb 
atoms is taken to be Qs — 102(6)ao [23], with oq the Bohr 
length; the characteristic length a±^ of the tight confinement 
in the y and z direction is = 57.6nm; the parameter e, in 
Eq.Q, characterizing the trapping in the axial direction, is 
given by £,■ = j dxVT{x)\^{x ~ Xi)\'^ ~ 8 x \Q-^ER{i~^Y, 
with ^'{x — Xi) the Wannier function centered around site / and 
L = 50 the total number of sites. The ratio U/J can be varied 
by changing the optical potential strength Vq. The temperature 
T and the number of particles in one tube are not directly ac- 
cessible experimentally. The averaging over an array of one- 
dimensional tubes in Ref. fl] can be understood as an averag- 
ing over condensates with different temperatures and particle 
numbers. However, one can understand the onset of the TG 
limit from simulations for a single tube with a fixed temper- 
ature. We used 7/7 = 1 at all interaction strengths, which is 
of the same order as the temperatures estimated in Ref. \^. 
In Fig. 12] the local densities and the momentum profiles are 
shown for several values of the optical potential strength Vb, in 
line with the actual values used in the experiment of Ref. 1^^. 
All Monte Carlo simulations consist of at least 20 chains of 
2'^ samples with each 50 — 200 off-diagonal updates such that 
error bars are not visible. 

Momentum profiles are experimentally measurable and can 
be calculated from a numerical simulation as 

n{p)^\^{p)\'Y.e-'P^J-'\b]b,), (4) 
jd 

where the envelope ^{p) is the Fourier transform of the Wan- 
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FIG. 2: Local densities n, in coordinate space as a function 
of the site index / (left) and the corresponding momentum pro- 
files n{p) as a function of the momentum p (in units hk) on the 
right. The axial optical lattice depths, the ratios U /J and the val- 
ues of the slope parameter a for soft-core (full line) and a' for 
hard-core bosons (dashed line) are (a) Vq/Ej^ = \,U /J = 1.75, a = 
2.71, a' = 1.69, (b) V^■,/ER = 5,(7/7 = 7.85, a = 1.92,a' = 1.38, 
(c) V^)|ER = 9.5, U/J = 28.6,a = l.OO.a' = 0.78, (d) Vq/Er = 
n,U/J = 52.28, a = 0.72, a' = 0.56, (e) Vq/Er = 20,U/J = 
258.54, a = 0.33, a' = 0.32. In each plot the average number of par- 
ticles is (A') ~ 15, the temperature isT/J = 1, and the lattice consists 
of 50 sites. 



nier function l'(x), p denotes momentum in units of hk and 
{b^jbi) is the one-body density matrix of the Bose-Hubbard 
model. In Fig. |5] the peak observed at p = 2hk is the first- 
order diffraction peak reflecting the presence of the optical lat- 
tice. The ratio between the height of the central peak and the 
first-order peak is solely related to the width of the Wannier 
orbitals and is not affected by averaging over the array of tubes 
or by the dynamics of the Bose-Hubbard model. The proce- 
dure to calculate the Wannier orbitals outlined above yields 
ratios in good agreement with the experimental data shown in 
Fig. 2 of Ref. |6]. This suggests that the ramping down along 
the axial direction in the experiment proceeded adiabatically, 



and it demonstrates that the discrete Bose-Hubbard model is 
a valid approach to describe the physics of ultracold atomic 
alkali gases in optical lattices. 

In every subplot (a-e) of Fig. |2l there is a region where 
the slope of the momentum distribution is almost linear (on 
log-log scale), similar to what occurs in an infinite homo- 
geneous Tonks gas at T = 0, which has an infrared diver- 
gence n{p) /?^'/^ at low momenta and an asymptotic tail 
n{p) o= p^'^ at high momenta 1 24]. In our case, the periodicity 
of the optical lattice sets an upper momentum scale pi^ = hJc. 
The width of the Wannier orbitals sets another upper scale, 
pw ~ (Vo/ERy^^pL, which turns out to be larger than pi for 
the parameter regimes considered here. The harmonic trap 
sets a lower momentum scale pr = (jnha))^^^ ~0.1/iA;, below 
which the momentum distribution is flattened because of the 
suppression of long-range correlations. Due to the trapping 
potential, the influence of temperature on the momentum dis- 
tribution will be different from the homogeneous case: ther- 
mal fluctuations will occur at the edges of the cloud and there- 
fore they will mainly affect the momentum distribution below 
Pt- Only the momentum distribution in the region between 
Pt and pi relates directly to the short-range dynamics of the 
Bose-Hubbard model and might show a power-law behavior 
similar to the homogeneous system. We have fitted the lin- 
ear parts of the log-log curves in this region with a power law 
n{p) o<: p^"". The slope a is sensitive to temperature, density 
and interaction strength, but to first order independent of the 
Wannier orbitals. 

By comparing the results for soft-core and hard-core bosons 
in Fig. |2j one sees that the TG regime is approached for op- 
tical potentials Vq/Er = 9.5 and Vq/Er = 12, while it is fully 
reached only at Vb /Er = 20, which in our model corresponds 
to a ratio U /J = 259. These values are in good agreement 
with Fig.^ where the energy for U/J — 200 is only 4% lower 
than the energy of an ideal Fermi gas. For lower optical po- 
tentials, finite boson-boson interactions certainly need to be 
taken into account and double occupancies in the center of 
the trap do play an important role. We see in subplot (e) of 
Fig. |2] that a Mott-like region is formed in the center of the 
trap. For a homogeneous system in the Mott phase, the dis- 
persion relation of the excitations has a gap of order U fTHl . 
meaning that the role of double occupancies is strongly sup- 
pressed 11211 . The Mott phase is entered at a ratio U/J as low 
as 1 .67 at r = for a density of one particle per site 1 25] . For 
the inhomogeneous system, the insulating behavior translates 
into a local compressibility that tends to vanish in the center 
of the trap \2M- Hence, for T ^U hard-core bosonic behav- 
ior can be reached for local densities varying from («,) = 
to (n,) = 1. However, the reduced local compressibility does 
not mean that for higher densities the TG regime would be 
reached at weaker optical potential strengths. Fig. |3l shows 
that a significant difference between the soft-core and hard- 
core momentum profiles persists even if the density-profile 
develops a Mott-like region, at an intermediate optical lat- 
tice strength Vq/Er = 7 {U/J = 14.3). This indicates that 
the short-range correlations in the Mott-like region differ sig- 
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FIG. 3: Density profiles and momentum density profiles n{p) 
(p in units hk) for different fillings A' at an optical lattice potential 
Vo/Er = 7 and temperature T/J = 1, for soft-core (full line) and 
hard-core bosons (dashed line). 



nificantly from the short-range correlations in the TG regime. 

In conclusion, we have shown that the experiment of 
Ref. 1 6] is very well described by a Bose-Hubbard model 
based on Wannier orbitals. Soft-core boson wave functions 
with a significant contribution of double occupancies can ex- 
plain the experimental results over the largest part of the opti- 
cal potential parameter range. Only for very deep optical lat- 
tices (Vq/Er = 20) do the atoms behave as hard-core bosons 
and does the Tonks-Girardeau picture apply. The averaging 
over the aiTay of one-dimensional tubes has only a minor ef- 
fect and does not significantly alter the momentum profiles. 
At higher densities, Mott-like correlations might develop, but 
they do not enhance the convergence to the Tonks-Girardeau 
regime. 
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